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We propose a new, self-consistent and dynamical scenario which gives rise to well-defined initial 
conditions for five-dimensional brane-world cosmologies with radion stabilization. At high energies, 
the five-dimensional effective theory is assumed to have a scale invariance so that it admits an 
expanding scaling solution as a future attractor. The system automatically approaches the scaling 
solution and, hence, the initial condition for the subsequent low-energy brane cosmology is set by the 
scaling solution. At low energies, the scale invariance is broken and a radion stabilization mechanism 
| drives the dynamics of the brane-world system. We present an exact, analytic scaling solution for a 

class of scale-invariant effective theories of five-dimensional brane-world models which includes the 
five-dimensional reduction of the Horava-Witten theory, and provide convincing evidence that the 
scaling solution is a future attractor. [HUTP-03/A066] 



o 



Q 

I. INTRODUCTION 

> 

Phenomenological models in which our universe is a four-dimensional (4d) hypersurface, or a 3-brane, embedded in 
' a higher-dimensional spacetime is currently of great interest. In the brane-world scenario, ordinary matter fields are 
t-H , confined to the brane, while the gravitational field can propagate in the extra dimensions (i.e., in the 'bulk') [1-5]. 
In particular, Randall and Sundrum proposed a self-consistent brane-world (RS2) scenario with 5d anti-de Sitter 
(AdS) bulk spacetime and showed that 4d gravity is also localized on the brane if the brane tension is positive and 
fine-tuned [5]. Many aspects of the Randall- Sundrum scenarios have been investigated, including weak gravity [6-8], 
the effective Einstein equation [9,10], cosmology [11-18], black holes [19-23] and so on. 

In the brane-world scenario, the evolution of our universe can be considered in two ways. On the one hand, it is the 
expansion of the induced geometry on our brane. Note that ordinary matter fields are assumed to be confined on the 
brane and, thus, propagate in the induced geometry. This picture is not very different from what we are accustomed to 
in the standard 4d cosmology. Indeed, the generalized Friedmann equation for homogeneous, isotropic brane universe 
in the RS2 scenario reduces to the standard Friedmann equation at low energy with corrections [11-15]. In this 
picture the bulk geometry also evolves 1 , corresponding to the expansion of our universe. On the other hand, the 
same evolution of our universe can be seen in a different way as a brane motion in a higher dimensional spacetime. In 
this picture the generalized Friedmann equation is understood as the equation of motion for the brane moving in the 
bulk spacetime. In the RS2 scenario this picture is remarkably simple for the homogeneous, isotropic brane universe. 
Indeed, the bulk spacetime is always locally AdS-Schwarzschild spacetime, irrespective of the brane motion [16-18]. 
In other words, the brane motion does not generate waves in the bulk. Unfortunately, in the more realistic scenarios 
which we shall consider in this paper this remarkable property is not satisfied. Nonetheless, it is still useful to keep in 
mind that there are two different ways to describe the same evolution of the brane-world cosmology: the expansion 
of the induced geometry and the brane motion in the bulk. 

Prior to the RS2 scenario, Randall- Sundrum had proposed a similar but different scenario called the RSI scenario [4]. 
This scenario involves two branes with fine-tuned tensions, one positive and another negative, and aims to solve the 
gauge hierarchy problem. In order to obtain an appropriate hierarchy between the Planck scale and the electroweak 
scale, the distance between the two branes must be set to about 50 times the bulk curvature scale. This itself is not a 
severe fine-tuning nor a problem at all, but it would be more satisfactory if this value could be dynamically realized. 
In the RSI scenario, unfortunately, from the 4d brane viewpoint the inter-brane distance represents a massless scalar 
degree of freedom and, thus, is arbitrary. Moreover, the existence of this massless scalar mode, called the radion, 
causes a more serious problem: the low energy 4d gravity on our brane is not general relativity but a Brans-Dickc 
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1 Here, the bulk geometry is supposed to be seen in the Gaussian normal coordinate based on the world-volume of our brane. 
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type theory with a too small Brans-Dicke parameter to be consistent with gravitational experiments [6]. Here, the 
radion plays the roll of the Brans-Dicke scalar. 

The problem of too large a deviation from the Einstein theory can be fixed if the radion obtains a non-zero mass. 
With a non-vanishing radion mass, the radion is not excited and becomes irrelevant at energies much lower than the 
radion mass, so that the 4d Einstein gravity is restored at low energy. This is called radion stabilization. Goldberger 
and Wise proposed a mechanism to stabilize the radion by introducing a scalar field in the bulk [24] . This scalar field 
has a potential in the bulk (bulk potential) and couples to branes via its potentials localized on the brancs (brane 
potentials). The introduction of scalar field(s) in the bulk is, of course, favorable from the viewpoint of string theory 
since a compactification from fO or 11 dimensions down to 5 dimensions in general introduces many 5d scalar fields. 
In this sense, the generalized RSI scenarios with radion stabilization seem more realistic than the RS2 model. With 
the Goldberger- Wise mechanism, 4d Einstein gravity is indeed shown to be recovered on our brane [25-29] 2 . The 
radion stabilization also helps recovery of the standard Friedmann equation at low energy [31,32]. 

The introduction of the bulk scalar complicates the brane-world scenario in (at least) two different ways. First, 
the gauge hierarchy problem is now entangled with the cosmological constant problem. The essential reason for this 
is that both the 4d cosmological constant and the stabilized value of the inter-brane distance depend on the bulk 
potential and the brane potentials of the scalar field in a non-trivial way. 

Second, the dynamics of the system becomes much more complicated. The bulk geometry is no longer simple if a 
brane is moving; i.e., if our 4d universe is expanding. A general brane motion generates waves of the bulk scalar field 
and makes the bulk geometry very complicated. This makes it extremely difficult to analyze the general dynamics 
of brane-world cosmology analytically. Recent development of BrancCodc [33] has made it possible to analyze the 
dynamics numerically, but there remains the problem of initial conditions. Namely, the system is so rich that it is 
not a priori trivial to choose physically relevant initial conditions for the numerical study. 

The purpose of this paper is to shed light on this second point in a particular way. In the rest of this paper, our 
main concern will be the dynamics of the brane-world system with radion stabilization. In particular, we shall propose 
a new dynamical scenario to give an initial condition for the brane-world cosmology with radion stabilization. 

II. PHYSICAL SCENARIO 

Suppose that the 5d theory has an approximate scale invariance at relatively high energy (but still much lower than 
the Planck scale) so that the 5d effective action allows a scaling solution. For example, we might imagine that the 
bulk potential V{4>) and the brane potentials \±{<j>±) have the form 

V(<f>) = V e- 2aK '^, A±(0±) = Aq [ e - aK *+± + e &K ^±] , (1) 

where <fi± is the pullback of 4> on each brane and a + (3± ^ 0. We assume that — (a + /3±)n 5 (j)± 3> 1 initially so that 
the term e^ ±K ^ ± can be neglected and the system has the scaling invariance. 

It is worth mentioning here that potentials of the form (1) arise in many theories of the fundamental interactions. 
For example, compactification of higher-dimensional theories, including superstring theories, gives rise to exponential 
terms in the effective potential. Terms with different exponents come from different physical effects, eg., higher- 
dimensional cosmological constant, curvature of compact manifold, anti-symmetric field flux, Casimir effects, etc. We 
shall provide more motivations for exponential potentials in subsection IV A. 

Since the inter-brane distance has the dimension of length, we expect that it should be proportional to time and, 
thus, expanding for the scaling solution. In this paper, we shall find such a scaling solution and present convincing 
evidence that it is a future attractor. For the moment, we just assume that the scaling solution is a future attractor. 
In this case, the system should approach it dynamically. Subsequently, since the Hubble expansion rate decays as 1/t 
and the size of the extra-dimension increases as t for the scaling solution, the system should enter a lower and lower 
energy regime. Hence, the scale invariance is no more than an approximate symmetry. In the above example (1), 
when — (a + f3±)n 5 (/)± becomes of order unity, the term e /3±/tB * ± in the brane potentials cannot be neglected any more 
and the approximate scaling invariance will be broken. 

Depending on the form of the effective action relevant to the lower energy regime, the inter-brane distance, called 
the radion, may play the role of an inflaton. Hence, we will have a 4d inflation on our brane driven by the radion. 
When the radion-induced inflation ends, the radion is stabilized and the fields confined on our brane can be reheated 
due to the oscillatory behavior of the system around the stabilized configuration. Subsequently, the conventional 



2 For an apparent conflict with the picture in ref. [30], and its resolution, see ref. [28]. 
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standard cosmology is realized. Since the radion is eventually stabilized, Einstein weak gravity is restored at low 
energy [25,27-29]. 

Hence, in this hypothetical scenario, the role of the scaling solution is to provide the initial condition for the radion- 
induced inflation and the subsequent evolution. To our knowledge there has been no well-defined scenario to provide 
an initial condition for the brane- world cosmology, except perhaps for the creation-from- nothing scenario [34] . If the 
scaling solution is a future attractor then the classical dynamics will automatically drive the system to the well-defined 
initial condition for the following evolution of the brane world cosmology. 

In the remainder of this paper we shall find an exact, analytic scaling solution and present convincing evidence for 
its attractor behavior. 



III. BASIC EQUATIONS 



In this section we present the basic equations for general potentials. In the following sections we apply these 
equations to a specific model with exponential potentials, motivated by string- inspired phenomenology 
In the 5d bulk we consider Einstein gravity and a scalar field: 
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d M <t>d M $-v($) 



(2) 



Hence, the 5d Einstein equation is 

Gmn = k 2 Tmn, 
Tmn — Om^nc, 



-d L <f>d L <b + v(4>) 



guN- 



The field equation for 



V 2 ^V'(0) = O, 



(3) 



(4) 



where V is the covariant derivative compatible with the metric gMN, follows automatically from the Einstein equation 
because of the Bianch identity V^G^ = 0. 

We assume that the extra dimension has the topology of Z2/S 1 and consider two end-of-the- world branes on the 
two fixed points. The tension A± of each brane in general depends on the pullback <j>± of the scalar field 4>. Hence, 
the brane action is of the form 



h = - dxt^qZX-^-) 



- j dx\ y / -q+ 



(5) 



where x± represents 4d coordinates on each brane and q± is the determinant of the induced metric. 

Note that we need to add the Gibbons-Hawking boundary term to the action, depending on the precise definition 
of I 5 . If J 5 includes the integration over the thin layers corresponding to the branes then the Gibbons-Hawking term 
must not be added since it appears automatically as we integrate the Einstein- Hilbert term over the thin layers [44]. 
On the other hand, if I 5 docs not include the integration over the thin layers then the Gibbons-Hawking boundary 
term must be added. In all cases, including the variation of the metric, the position of the hypersurface and the scalar 
field, the variational principle gives the correct set of equations of motion: the Einstein equation and the field equation 
of the scalar field off the branes, Israel's junction condition and the matching condition for the scalar field [44]. 



A. Bulk equation 

For simplicity we assume that the geometry on each brane is everywhere and for all times described by the 4d flat 
FRW universe and that the two branes are parallel to each other. In this case, it is natural to expect that the 5d bulk 
spacetime between the two branes also has the same symmetry. Hence, we assume that the metric and the scalar field 
have the following forms: 

ds 2 = -n(t,y) 2 dt 2 + a(t,y) 2 5 ij dx i dx j + b(t,y) 2 dy 2 , 

<t> = <t>{t,y), (6) 
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where i = 1,2,3; y represents the extra dimension; and the world- volume of each brane is expressed as 

y = Y ± (t), (y_(t)<y+(t)). 



For this ansatz, we obtain 
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where a dot and a prime denote derivatives with respect to t and y, respectively. 
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B. Boundary condition 



On each brane at y = Y±(t), the induced metric q±p V , the extrinsic curvature K± jlv and the pullback <f>± of the 
scalar field play important roles. They are given by 



q±u 



n 2 - b 2 Y? 



q±ij = a 2 Sij, 
K± = Kl t = 



bY + 
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Having the Z 2 symmetry, the Israel junction condition and the scalar field matching condition are written as 
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£± = ±-g-A±(^±), 



£±=±yA±(0 ± ), 
1 



When Y± (t) are constants y± , these conditions reduce to the following simple conditions 
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IV. SCALING SOLUTION 



In this section we shall seek a scaling solution for a specific model with exponential potentials. For a special choice 
of parameters, the model reduces to the 5d reduction of the Horava-Witten theory [45]. On the other hand, the 
solution we shall find does not appear to reduce to the known flat FRW solution in the Horava-Witten theory [46]. 
(For closed and open FRW solutions in Horava-Witten theory, see rcf. [47].) 



A. Action 



Exponential potentials of the form 



V = V e 



-2and 



(14) 



for (dilaton coupling) constant a, arise in many theories of the fundamental interactions including superstring and 
higher-dimensional theories. Typically, 'realistic' supergravity theories predict steep exponential potentials [37] (i.e., 
2a 2 > 1). The effective action of the Horava-Witten theory has a 2 = 2 (see below). A smaller effective a can 
arise in assisted theories which contain more than one scalar field [36]. In addition, for many compactifications 
of higher-dimensional theories there is a consistent truncation to a single scalar field <j> an d positive potentials of 
this form which arises via generation by non-zero flux of antisymmetric tensor fields (the T 7 compactification of 
lid supergravity with non- vanishing 4- form field strength yields an exponential potential with a = \/7; in general, 
flux compactifications seem to yield a > \/3 [38]), and hyperbolic compactifications in the context of String/M- 
theory [39,40] (the compactification of lid supergravity on a 7d compact hyperbolic space has a = 3/\/7; in general, 
hyperbolic compactifications seem to lead to 1 < a < y/3 [41]). A comprehensive qualitative analysis of scalar field 
cosmological models with an exponential potential has been presented [42]. 

In this paper we consider a brane-world analogue of such models with exponential potentials. To motivate our model 
with exponential potentials, let us briefly review the 5d reduction of the Horava-Witten theory. The 5d effective action 
for the Horava-Witten theory is [45] 
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(15) 



where $± is the pullback of the bulk scalar field $ onto each brane. If we introduce a canonically normalized scalar 
field <f> by 



V~2K 5 (j) = In 



(16) 



this action is reduces to 
Ihw = / dx 5 v^g 



-2^/2K 5 q 



>M<1 
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where <j>± is the pullback of the scalar field <j> on each brane. 
The model we shall consider is 



dx b 



-J dxt^TX^e-^ 



J dx 4 + ^^X+e- aK ^, (18) 



where a, V and are constants. This action includes the 5d effective action of the Horava-Witten theory as a 
special case (a = \/2)- 



B. Ansatz 



We would like to seek a special solution which corresponds to an 'equilibrium state', which represents a future 
attractor for the model. In a general situation where branes move arbitrarily, the motion of branes produces waves 
of the scalar field <f> in the bulk. The waves in the bulk interact with both branes in the sense that they can come 
and go between two end-of-the-world branes. Hence, the general situation is very complicated, and non-local from 
the viewpoint of a 4d observer on a brane [48] . On the other hand, we would like to seek a special solution in which 
waves emitted from one brane and those from another brane are in equilibrium so that there is effectively no scalar 
wave in the bulk. Hence, what we would like to seek is a very special situation in which the non-local effects due to 
scalar waves are completely suppressed. 

In this kind of 'equilibrium' situation without waves in the bulk, nothing should propagate from one brane to another 
and, thus, we expect that all physically meaningful functions of t and y should be a product of a function of t and a 
function of y. In particular, in this case the time dependence of n can be removed by a coordinate choice. Moreover, 
from our experience in 4d scalar field cosmology with an exponential potential, we expect that the 'equilibrium' 
situation should correspond to a power-law expansion on each brane, a(t,y±) oc t p [42]. Actually, in 4d cosmology 
with an exponential potential, a power-low expansion is an attractor of the system. One of the essential physical 
reasons for this is that all physically relevant, dimensionful quantities scale in tandem, according to dimensionality. 
For example, H 2 oc k 2 V oc n\pki n oc 1/t 2 , where H is the Hubble expansion rate, K4 is the 4d gravitational coupling, 
V is the exponential potential and pkin is the kinetic energy of the scalar field. In our model, besides the Hubble 
expansion rate on each brane, the inter-brane distance is another physically relevant, dimensionful quantity. Since it 
has the dimension of length, we expect it should also be proportional to t in the 'equilibrium' situation. Thus, not 
only the 4d universes on branes but also the inter-brane distance should be expanding. 

For the reasons explained above, we consider the following ansatz. 

ds 2 = n{y) 2 (-dt 2 + t 2p 6 ij dx i dx j + t 2 dy 2 ) , 
p -2</>(y) 

K 5 \v \e - p f2 , 

Y±(t)=y±, (19) 

where n(y) and i[>(y) are functions of y, and p, I and y± (y_ < y + ) are constants. Note that t, y and x l are 
dimcnsionless and that n{y) and I have the dimension of length. 

C. Equations 

With this ansatz, the Einstein equation in the bulk reduces to 

3a 2 p -1 = 0, 
(]nn-pipy = 0, 

vt-^vtrsF''-*- - (20) 

where the plus and minus signs correspond to Vq > and Vq < 0, respectively. Hence, we obtain 



n = noe 1 " 1 ', (21) 



G 



where no is an arbitrary constant and i\> is a solution to the following equation. 

(^') 2 - 1 ± e 2(p_1) * = 0, (22) 

where we have set 



'=yj43&?r a «o. (23) 

and the plus and minus signs correspond to (4 — 3a 2 )Vb > and (4 — 3a 2 )Vo < 0, respectively, and we have assumed 
that a 2 ^ 4/3. 

The Israel junction condition (12) is written as 



e -(p-i)iy 



c^KsApno _ /|4-3a 2 | 



= ±7±, 7± = T^- = ——j==\[ 1 " 1 . (24) 



The scalar field matching condition (13) is actually the same as the above junction condition and, thus, does not 
provide any independent boundary conditions. 

The physically relevant quantities are the power index p = of the expansion of the branes, the warp factor 
W = n(y+)/n(y-), and the ratio v± of the inter-brane distance L = t n{y)dy to the proper time r± = n(y±)t on 
each brane at y = y± . These are given by 

W = n ^+} = I± = e v[Hy+)-^{y-)\ v± = J± = e -piKi/±) f V+ e P^(v)dy (25) 
n(y-) r- t± J v _ 



D. Solutions 

For (4 — 3a 2 )Vb < 0, the solution to the bulk equation is 

(p-i)v = / sinh [\p -M(y- vo)] ( for yo < y) 

sinh[|p- l\(y - y)] (for y < y Q ) ' 



(26) 



where yo is a constant. Note that tp' diverges at y = yo. Hence, yo must not be between the two branes. For this 
solution, the junction condition reduces to 

cosh[(p - l)(y± - yo)} = T7± (for yo < y- < y+) ^ 
cosh[(p- l)(y± -j/ )j = ±7± (for V- < V+ < 2/0 ) 

Hence, if and only if 1 < 7_ < —7+ or 1 < 7+ < —7-, the junction condition uniquely determines y± — yo- The 
physically relevant quantities W and v± are given by 

W=(^±Y^, , ± = (7|-D^r * (28) 

\r+- l J h- [smh\(p-l)(y-y )\}—i 

Not only W but also v± are uniquely determined by p (or equivalently \a\) and 7±. For example, when p = 2, 

V^5T= 17+1^71 - 1 + 17-1^7+ - 1- (29) 
For (4 — 3a 2 )Vo > 0, the solution to the bulk equation is 

e -CP-i)* = cosh [(p- l)(»-i/ )], (30) 
where yo is a constant. For this solution, the junction condition reduces to 

sinh[(p -l)(y±-y )} = T7± • (31) 
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Hence, if and only if 7_ + 7+ < then the junction condition uniquely determines y± — yo. The physically relevant 
quantities W and v± are given by 

W v ± = {l l + l)^r * r . (32) 

W+ + l J Jy- [cosh|(p-l)(y-y )|]~ 

Again, not only W but also v± are uniquely determined by p (or equivalently \a\) and -f±. For example, when p = 2, 



V^T^=-7+V7- + l-7-V^+ + L ( 33 ) 



V. STABILITY AGAINST LINEAR PERTURBATIONS 



Now let us argue that the expanding scaling solution is a future attractor. For this purpose we investigate linear 
perturbations around the scaling solution and show stability. The linear perturbations analyzed in this section include 
all important physical effects such as boundary conditions on the branes, fluctuations of brane positions and scalar 
waves in the bulk. The bulk scalar waves generated by fluctuations of brane positions threaten to destabilize the 
expanding scaling solution, but we shall explicitly sec that this is not the case and that the scaling solution is stable. 
In the following, for simplicity we consider the case where ip'(y) 7^ for y_ < y < y + . 

Let us consider linear perturbations around the scaling solution found in the previous section. 

ds 2 — —n(t, y) 2 dt 2 + a(t, y) 2 5ijdx l dx : > + b(t, y) 2 dy 2 + 2e(h ty dtdy + hudtdx 1 + h yi dydx l ), 
n(t,y)=n e^M [1 + e 5n(t,y)} , 
a(t,y)=t"no^ [1 + eSa(t, y)} , 
b(t,y)=tn e p ^ [1 + e6b(t, y)] , 

an b ^{t, y) = In [k 5 /T^ **] + My) + ety>(t, y), (34) 
where p and if>(y) are given in the previous section. From the symmetry, we set 

h ty = h ti = h yi = 0. (35) 

By introducing a new time coordinate 

r = lnt, (-oo<t<oo), (36) 

it is shown that all coefficients in all relevant equations (i.e., the Einstein equation, the Israel junction condition and 
the scalar field matching condition) linearized with respect to e are independent of r. Thus, it is convenient to Fourier 
expand the perturbations as 



5n(t, y) = 5n(y)e 
5a(t,y) = Sa(y)e~ 
5b(t,y) = 6b(y)e~ 

6iP(t,y) = 6j(y)e- i ^. (37) 



-2UJT 

? 

-iur 



A. Gauge conditions 

The set of variables (h tyi u,yi, Sn, 5a, Sb, 8ip) includes not only physical degrees of freedom but also gauge degrees 
of freedom. The infinitesimal gauge transformation is 

50 ^ 50 - e^V (38) 
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where V is the covariant derivative compatible with the background metric, <^°) is the background of <f> and is an 
infinitesimal vector representing the gauge degrees of freedom. Hence, each component transforms as 

8n->5n+\ [(1 - iw%{y) - p^'{y%{y)] , 
6a -> 5a + \ [p|t(y) - pip' \y%{y) - £\\{y)} , 



n 



o 



h ty ^h ty + te 2 ^ [iwUy)-Uy)} , 
h ti - &« + iort 2 "-^^^^ (y)<W - 4), 

^ - h yi - t 2 *e 2 ^% (y)5 ij (x j - 4), 

6j^Sj + ±[£ t (y) - i>'{y%{y% (39) 
where we have Fourier expanded £^ as 

i v =t 2 e 2 ^y%{y)e-^, 

e . = t2p e 2p*(w)^ { y )e-^Sij{aP - 4). (40) 

Here, Xq (j = 1,2,3) are constants. The symmetry assumption h ty = h ti = h y i = is consistent with the gauge 
transformation if and only if 

i^ v (y) - Uv) = »k (y) = £] (w) = o. (41) 

By using the gauge freedom, we can set 

Stp = 0. (42) 

However, this condition does not fix the gauge freedom completely. So, let us investigate the residual gauge freedom. 

(i) For w^O and — ioj ^ p — 1, the residual gauge freedom is 

Sn —> Sn - IUJ+ P — -£(?/), 
6a — > <5a, 

«-«-^ <«) 

where £ (= £ t ) is an arbitrary solution of 

^'f = (44) 

Here, we have used the background equation (22) and its derivative. 

(ii) For w^O and — iuu = p — 1, there is no residual gauge freedom. 

(iii) For u> = and p ^ 1, the residual gauge freedom is 

Sh — > <5n + <5rio, 
<5a — > <5a + oao, 

(45) 

where <5no and Sao are constants. 

(iv) For u = and p = 1 , the residual gauge freedom is 

5a — > (5a + <5ao, (46) 

where (5ao is a constant. 
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B. Linearized equations 



By using r as a time coordinate, the Einstein equation linearized with respect to e becomes relatively simple in 
the sense that there is no explicit r dependence in any coefficients. The Israel junction condition and the scalar field 
matching condition also share this property. This is the reason why we have Fourier-expanded the linear quantities 
(5n 5a, 5b, 5ip) with respect to r. Hence, what we have is just an eigenvalue problem in one-dimension. 

The (Id) bulk equations are 

^ s _, = _ uj 2 + (p- ljiw s _ + (p _ ( ^)2 ^ + ^ _ ^ _ ^ 
ip'Sa' — —iujSa + p {ip') 2 5b — p5n, 

p L J 

where the plus and minus signs correspond to (4p — l)Vo > and (4p — l)Vo < 0, respectively. 

Since the gauge freedom has been fixed only up to the residual gauge freedom, the brane position is not at y = y± 
any more. Hence, we need to consider the perturbed position of the brane: 

y = Y ± (r) = y ± +e5y ± e-^ T . (48) 

For the fluctuating position of the brane, the formulas (12) and (13) cannot be used. As explained in Appendix A, 
the general boundary condition (11) gives the following boundary condition for perturbations: 

5h' — iuoxp'Sb, 
5a = 0, 

(p - 1 + iu)6y± = -ip'Sb. (49) 
Note that these boundary conditions are imposed on y = y± (not Y±(y)). 



1. Second-order equation for 5a 

From the set of equations (47) and the background equation (22), we can show that 5a satisfies the following 
second-order ordinary differential equation: 

5a" + Spip'Sa' + lu(lo + 3ip)5a = 0. (50) 

Now let us show that 3o> < unless 5a = 0. By using the equation (50) with the boundary condition 5a! — at 
V = V±i h is easy to show that 

Lu(uj + -iip)l dye 3p ^\5a\ 2 = dye^\5a'\ 2 . (51) 
Jy- Jy- 

Hence, ui(u) + 3ip) is real and non-negative unless 5d = 0. The eigenvalue lo(uj + 3ip) can vanish only if 5a! = 0. Thus, 

-3p < 3?w < 0, (52) 

or 

uj = 0, 5a = 0. (53) 

In the latter case, 5d can be set to zero by using the residual gauge freedom (45) or (46). Therefore, we have shown 
that -3p < 3?w < or 5a, = 0. 
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2. Solution with 8a — 

When 5a = 0, the set of first-order equations (47) reduces to 

ip 5h = iujSn, 

{^') 2 5b = 5n. (54) 

For — iuj p — 1, the boundary condition (49) is automatically satisfied, but any solutions to (54) with 5a = is pure 
gauge because of the residual gauge freedom (43) or (45). On the other hand, for — ioj = p — 1, there is no residual 
gauge freedom, but the third boundary condition in (49) becomes 5b — at y = y± and, thus, (54) implies that 
5b = 5h = 0. Therefore, we have shown that there is no non-vanishing physical solution with 5a — 0. 



3. Stability 

In summary, we have found that there is no unstable mode with 5a ^ and that there is no physical mode with 
5a = 0. Thus, the scaling solution is stable against linear perturbations. 



VI. TOWARDS A NON-LINEAR STABILITY ANALYSIS 

Following the linear perturbation stability analysis, we need to establish the stability at the non-linear level. In 
general this is very complicated. However, we can investigate some of the non-linear effects by considering a subclass of 
models with a particular ansatz. This will give us an indication of the qualitative features of the non-linear dynamics. 
We shall perform a more comprehensive analysis at a later time. 

The ansatz we consider is 



ds 2 = n 2 Q e 2p ^ [-N(t) 2 dt 2 + A(t) 2 5 l3 dx l dxi + B(t) 2 dy 2 ] , 



(55) 



where n is a constant with the dimension of length. The scaling solution is, of course, included within this ansatz. 
The (ty)-component of Einstein equation leads to 



0(t) = In[B(t)] + Go, 



(56) 



where 0o is an arbitrary constant. Since <do can be absorbed by a redefinition of ip(y), for convenience we can choose 
0o as 



whence the field equation reduces to 



O 26 = ^4n 2 \Vo\ 
3p|4p-l|' 



(57) 



(l + 2p)^-2(p-l)^- 
(l + 2p)| + 3p|J 



(V0 2 ±e 2{ P-^- 



B N 
(4p-l)- + (2p+l)- 



A N 
3(2p-l)--(2p+l)- 



A 2 B 2 
B A 2 

b- & T 2 = {) ' 



p{2p+l)N 2 



A 2 B 2 
' A 2 



ABB 
2 —- pB 



0, 



(58) 



where the plus and minus signs in the last equation are for (4p — l)Vo > and [Ap — l)Vb < 0, respectively. Hence, 
the model is separable. 

Hereafter, we set N(t) to a constant by choice of time coordinate. Defining 



A h B 



(59) 
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the last equation in (58) becomes 

clB- 2 = 2a 2 + 2ab - pb 2 , (60) 

where c 2 , is effectively the rescaled separation constant, and the evolution equations are then 

(1 + 2p)a = -3a 2 + {Ap - l)ab - 3p 2 b 2 , (61) 
(1 + 2p)b = 6a 2 - 3{2p - l)ab - (1 + bp)b 2 . (62) 

Differentiating equation (60) and using equations (61) and (62) we get an expression which is satisfied identically, so 
that (60) is a constraint (that propagates along the solution curves). 

We first note that since the system (61), (62) is homogeneous, we can define 

x = a/b, b = Inb, (63) 

and the system reduces to a single ordinary differential equation: 



dx _ 6(x -p)(x + i(l + ^/T+2p))(x + \{1 - VTT2^)) 
db ~ (1 + 5p) + 3(2p - l)x - 6x 2 

which can be integrated to obtain 



(64) 



b = b 2 (x-pf(x + i(l + VT + WOr + \{l - JT+2~p))\ (65) 

where 6g is an integration constant and 
1 

a = — , 
3p 

~ -3 P - 1 + y^TT (3V2^TT+l)(V2^TT-i) = (3V27TT + 5)(V2^TT-i) 

/ ~ 6p 12p 12p 



. = -3p - 1 - V2p-+T = (3727+T-l)(V2^+T+l) = _ (3V2j7+T - 5) (VW+T + 1) 
7 6p 12p 12p ' 

Note that 

a + /3 + 7 = -1, 
d > 0, 
-1< £ < 0, 

7 < (67) 

for any positive p = l/(3a 2 ), and 7 > —1 if p > 8/9. 

We can analyze the asymptotic form of (65), but perhaps a better way to show the late-time stability of the scaling 
solution is as follows. We rewrite equation (60) as 



(2p+l)b 2 egg-* 



Defining 



Z ~ 2a + b ~ 1 + 2x ' lbyj 



we see that z is bounded; 1 — z 2 > 0. Using (65), the evolution equations (61) and (62) then become a single ordinary 
differential equation for z: 

36 ° : (l-yiT2^r +1 (l + z) /5+1 (l-zF +1 . (70) 



2 2 ' 3 V1 + 2p 
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Since z 2 < 1, this constitutes a one-dimensional non-linear dynamical system. 
The exponent (a + 1) is positive definite, so that 



* - 7rm (71) 

(where < z s < 1), is always an equilibrium point of (70). Depending on the signs of the exponents (,9 + 1) and 
(7+1), z± = ±1 are also equilibrium points (z_ is always an equilibrium point, while z + is an equilibrium point for 
p > 8/9). 

The solution of (70) close to z s is given by 

z = z s -z 2 (t + t )- 1 / & , (72) 

so that z — > z s as t — > 00 since a > 0. Hence, z s is a local sink. Since z is bounded, z s is a global attractor in the 
physical phase space. Indeed, since {3 < and 7 < 0, z± act as local sources and all physical solutions evolve from 
one of z + and z_ to z s . 

From equations (69) and (71), the global attractor has 

- b {=x)=p (73) 



so that equation (60) yields 



and hence after a time translation 



B 2 = r / a - = const. (74) 
p(2p+l) V ' 



B = B t, A = A t p , (75) 

so that the global attractor is the scaling solution. 

To connect with earlier work, the approximate solution of (61) and (62) at late times (i.e., the linearized solution 
around the attractor (75)) is (N = 1): 



A = A t p [l + Zpdt- 1 +c 2 t- 3p 

Cl ,_i c 2 



B = B t 



1 _ ^t-i _ ^ t - 3p 
p p 



(76) 



where the Cj are arbitrary constants (and the constants are subject to the constraint (60); e.g., Cq = p(l + 2p)B 2 ). 
From this we immediately sec the decaying modes and the local stability of the attractor. 

Consequently, we have an approximate bulk (5d) solution at late times. The 5d Ricci tensor is given by its attractor 
values of order 0(i -1 ) and 0(t~ 3p ) (leading to pk in <~ t~ 2 ), and the 5d conformal tensor (which is zero for the exact 
attractor solution) can be evaluated to leading order. We can therefore calculate the irreducible decompositions U, 
Q^, and of the projected bulk Weyl tensor £ M „ [9,10], from which we find that 

Qn = 0, Vfa, = (77) 

and 

U ~ [0(t- 3p ) + 0{t- 2 )} Pkm « Pkm . (78) 



VII. DISCUSSION 



We have investigated a class of scale-invariant effective theories of 5d brane-world cosmology with a bulk scalar 
between two end-of-the- world branes. As a special case, this class includes the 5d reduction of Horava-Witten theory. 
We have found an exact, analytic scaling solution in which the scale factor of our 4d universe and the inter-branc 
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distance expand as t p {p > 0) and t, respectively 3 . It is perhaps worth mentioning that the scaling solution with 
p > 1 corresponds to a power-law inflation on our branc. The scaling solution corresponds to an equilibrium state in 
which the motion of the branes does not produce any scalar waves in the bulk. Because of this remarkable physical 
property the scaling solution is expected to be a future attractor of the system. Indeed, we have presented convincing 
evidence for the attractor behavior: the stability of the scaling solution against general linear perturbations with the 
4d FRW symmetry and stability within a class of non-linear perturbations. 

Based on the attractor behavior of the scaling solution, in Sec. II we proposed a new, self-consistent and dynamical 
scenario in the early brane-world universe. First, since the scaling solution is a future attractor, the system is 
automatically driven towards it as far as the effective action at high energy is of the scale invariant form. Second, as 
the energy scale becomes sufficiently low according to the evolution along the scaling solution, the scale invariance 
of the effective action should be broken at some point. After that, the brane-world system deviates from the scaling 
behavior and starts to be driven by the radion stabilization. The radion stabilization guarantees that 4d gravity on 
our brane is described by Einstein theory and the standard cosmology is realized at low energy. 

In this scenario, the scaling solution plays a central role: the attractor behavior of the scaling solution makes it 
possible to give well-defined initial conditions both in the bulk and on the brane for the evolution after the breaking 
of scale invariance. It is also the scaling solution that brings the system from high energy to low energy both in the 
bulk and on the brane and, thus, triggers the breaking of the scaling invariance. 

Having a well-defined initial condition given by the scaling solution, a natural question arises: "What kind of 
evolution can we expect subsequently?" One interesting possibility is an epoch of inflation driven by the radion (i.e., 
the inter-brane distance). As stated above, after the scale invariance of the effective theory is broken at some point, 
the radion stabilization mechanism takes over and starts driving the brane-world system. Indeed, we have suggested 
a simple illustrative action in which the transition from the scale-invariance to the radion stabilization is smooth. 
During the transition, the dynamics of the bulk scalar field and the radion can drive the evolution of our 4d universe. 
Of course, even without the scaling behavior at high energy, the bulk scalar and the radion can drive the system. 
A big difference, however, is that one cannot expect a well-defined, predictable initial condition in the case without 
the high-energy scaling behavior since the evolution from an arbitrary state to the radion stabilization should be a 
very violent process and should involve fully non-linear (scalar and gravitational) waves in the bulk. Our scenario 
based on the high-energy scaling behavior makes the evolution towards the radion stabilization much smoother and, 
more importantly, predictable. In particular, in our scenario, there may be an inflation in the transition epoch. This 
possibility is certainly an interesting subject for the future work. 

Another question that might naturally be asked is: "What is the beginning of the brane-world before the scaling 
behavior?" Actually, since the attractor behavior of the scaling solution would make the low-energy evolution of the 
brane cosmology almost insensitive to the beginning, this question might be irrelevant for our scenario. Nonetheless, 
it may be still interesting to think about this kind of question. As already explained, for the scaling solution both 
the scale factor of our universe and the inter-brane distance expand. This implies that both our 4d universe and the 
extra dimension were extremely small at an early epoch. This is a completely trivial statement in our scenario, but it 
would not be so trivial if we did not have the high-energy scaling behavior. Without the high-energy scaling behavior, 
the initial value of the inter-brane distance for the evolution towards the stable value can be either small or large. On 
the other hand, in our scenario the initial value of the inter-brane distance should start from a small value. Is this 
favorable from the viewpoint of quantum gravity? Can we consider the smallness of the initial inter-brane distance 
as an indication of a brane collision or brane scattering? 4 The last two questions are evidently outside the scope of 
this paper but will be interesting future projects. 

It is perhaps worthwhile asking whether we can expect similar scenarios to work in more general situations and/or 
theories. When the dynamical effects of fields are negligible (eg., the bulk waves), then we might expect the same 
qualitative behavior and hence the scenario to work. For example, we expect predictable, smooth evolution towards 
the radion stabilization to be valid in more general situations whenever a symmetry at high energy exists, leading to 
attractor behavior, which is subsequently broken by radion stabilization at low energy. 

As already stressed many times, we have presented convincing evidence that the scaling solution is a future attractor. 
Since the attractor behavior is so essential in our scenario, in future work we shall further investigate the full non-linear 
dynamics both analytically and numerically. We shall also include other matter fields on the branc. Since the scale 
invariant theory is motivated by compactification of higher dimensional theories (see Sec. II and subsection IV A), it 



3 After submission of this paper, the authors were informed of ref. [49], in which a closely related solution had been found. 
4 It is perhaps interesting to note that the expanding scaling solution we found is a conformally Kasner geometry and that 
Kasner-like geometries were found to be generic collapsing solutions in brane collisions [33] . 
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is also worth investigating highcr-dimcnsional interpretations of the scaling solution. In particular, it is interesting to 
generalize the scaling behavior to higher dimension (greater than 5) and study the resulting physical consequences. 
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APPENDIX A: BOUNDARY CONDITION FOR LINEAR PERTURBATIONS 

In this appendix we explain the derivation of the boundary condition (49) for linear perturbations. For more general 
prescription, see ref. [50]. For general trajectories of branes the boundary condition is given by (11) at y = Y±(t). For 

n(t, y) = n e p ^ [l + eSn(y)e- WT ] , 
o(t, y) = PnoeT+W [l + e5a(y)e-^ T ] , 
b(t, y) = trioe 13 ^ [l + eSb(y) e - luJT ] , 
aK 5 <t>(t,y) = ln k 5 */\Vo\ It + 4>{y), 

Y±(T)=y± + e5y±e- i " r , (Al) 
the extrinsic curvature components JC and K,, and the normal derivative d±<p± of the scalar field are expanded as 

K± =K,f +e/ci 1} +0(e 2 ), 
K± = £i 0) + e^£ 1) e- <WT + 0(e 2 ), 
d±.4>± = d^f + ed ±( f> { f>e- iuT + 0(e 2 ), (A2) 



where 



n t 



d±4>± = t, A3 

otK^not 



and 



n t 

JC + = I"— ipu)8y± + 8a! — pip' 5b] 



n t 

p-pip 

d±4>± = 7 [-iuSy ± - iP'Sb] . (A4) 

an^not 

What is important here is that the boundary condition in the form (11) must be imposed on the perturbed position 
y = Y±(t) of the branc. On the other hand, it is mathematically convenient to impose a boundary condition on a 
fixed position in the coordinate y. Hence, let us convert the boundary condition (11) at y = Y±(t) to a boundary 
condition at y = y± . The result up to the linear order in e is 

IC ( i Y Sy ± + = ±^(d 4)± \ ± )ip , 5y±, 
(9±4 0) )'%± + &±4> i ± ) = T^{dl ± \±W8y ± . (A5) 
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This form of the boundary condition must be imposed at y = y± . Substituting the above expressions for i^, 

and d±</J± and simplifying the expressions by using the background equations, we obtain the boundary condition 
(49). 
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